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Abstract 

This paper addresses the set-point control problem of a heat equation with in-domain 
actuation. The proposed scheme is based on the framework of zero dynamics inverse 
combined with flat system control. Moreover, the set-point control is cast into a mo¬ 
tion planing problem of a multiple-input, multiple-out system, which is solved by a 
Green’s function-based reference trajectory decomposition. The validity of the pro¬ 
posed method is assessed through the convergence and solvability analysis of the con¬ 
trol algorithm. The performance of the developed control scheme and the viability 
of the proposed approach are confirmed by numerical simulation of a representative 
system. 

Keywords: Distributed parameter systems; heat equation; zero-dynamics inverse; 
differential flatness. 


1. Introduction 

Control of parabolic partial differential equations (PDEs) is a long-standing prob¬ 
lem in PDE control theory and practice. There exists a very rich literature devoted to 
this topic and it is continuing to draw a great attention for both theoretical studies and 
5 practical applications. In the existing literature, the majority of work is dedicated to 
boundary control, which may be represented as a standard Cauchy problem to which 
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functional analytic setting based on semigroup and other related tools can be applied 
(see, e.g., mm®®). It is interesting to note that in recent years, some methods that 
were originally developed for the control of finite-dimensional systems have been suc- 
10 cessfully extended to the control of parabolic PDEs, such as backstepping (see, e.g., 
GOB El), flat systems (see, e.g., lf8l l9l ITOl ITTl 12, 13, 1_4|), as well as their variations 
(see, e.g., HgJITB, 171 ). 

This paper addresses the in-domain (or interior point) control problem of a heat 
equation, which may arise in application related concerns for, e.g., the enhancement of 
is control efficiency. Integrating a number of control inputs acting in the domain will lead 
to non-standard inhomogeneous PDEs HUGH, which should be treated differently than 
the standard boundary control problem. The control scheme developed in the present 
work is based-on the framework of zero-dynamics inverse (ZDI), which was introduced 
by Byrnes and his collaborators in lH9l and has been exploited and developed in a series 
20 of work 120 21 22]£23i'24|. It is pointed out in |[23l that “for certain boundary control 
systems it is very easy to model the system’s zero dynamics, which, in turn, provides 
a simple systematic methodology for solving certain problems of output regulation.” 
Indeed, the construction of zero dynamics for output regulation of certain in-domain 
controlled PDEs is also straightforward (see, e.g., j22l ) and hence, the control design 
25 can be carried out in a systematic manner. A main issue related to the ZDI design is that 
the computation of dynamic control laws requires resolving the corresponding zero dy¬ 
namics, which may be very difficult for generic regulation problems, such as set-point 
control. To overcome this difficulty, we leverage one of the fundamental properties of 
flat systems, that is if a lumped or distributed parameter system is differentially flat (or 
30 flat for short), then its states and inputs can be explicitly expressed by the so-called flat 
output and its time-derivatives mm. In the context of ZDI design, the control can 
be derived from the flat output without explicitly solving the zero dynamics. Moreover, 
in the framework of flat systems, set-point control can be cast into a problem of motion 
planning, which can also be carried out in a systematic manner. 

35 The system model used in this work is taken from |221 . In order to perform control 
design, we present the original system in a form of parallel connection. This formu¬ 
lation allows a significant simplification of computation. As the control with multiple 
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actuators located in the domain leads to a multiple-input, multiple-output (MIMO) 
problem, the design of reference trajectories is not trivial. To overcome this problem, 

40 we introduce a Green’s function-based reference trajectory decomposition scheme that 
enables a simple and computational tractable implementation of the proposed control 
algorithm. 

The remainder of the paper is organized as follows. Section[2]describes the model 
of the considered system and its equivalent settings. Section [3] presents the detailed 
45 control design. Section [4] deals with motion planning and addresses the convergence 
and the solvability of the proposed control scheme. A simulation study is carried out 
in Section[5] and, finally, some concluding remarks are presented in Section[ 6 ] 

2. Problem Setting 

In the present work, we consider a scaler parabolic equation describing one-dimensional 
heat transfer with boundary and in-domain control, which is studied in ll22l . Denote 
by z(x, t) the heat distribution over the one-dimensional space, x, and the time, t. The 
derivatives of z(x, t ) with respect to its variables are denoted by z x and z t , respectively. 

For notational simplicity, we may not show all the variables of functions if there is no 
ambiguity, e.g., z = z(x, t). Consider m points Xj, j = 1 , ..., m, in the interval ( 0 , 1 ) 
and assume, without loss of generality, that 0 = x$ < X\ < x^ < ■ ■ ■ < x m < 

m 

,T m+ i = 1 . Let f2 = (J (xj,Xj+ 1 ). The considered heat equation with boundary and 

i =o 

in-domain control in a normalized coordinate is of the form: 


Zt(x, t ) - z xx (x, t) = 0, x G tt, t > 0, 

(la) 

z(x,0) = <j>(x), 

(lb) 

o' 

II 

HO 

cT 

o 

1 

Ho' 

o' 

II 

o 

0) 

(lc) 

Biz = z x { 1, t) + hz{l, t) = 0, 

(Id) 

z(xj) = z(x~), j = 1,..., m, 

(le) 

B X j z = = Uj (t ) , j 1,2,..., 771, 

(If) 
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where for a function ip(-) and a point x £ [0,1] we define 

bP)x = tp(x + )-tp(x~), 

with x~ and x + denoting, respectively, the usual meaning of left and right hand limits 
to x. The initial condition is specified by ([Tb]) with 4>{x) £ L 2 ( 0,1). In System <[l), we 
assume that we can control the heat flow in and out of the system at the points Xj, i.e., 

Uj(t) = [z x ] x . = z x (xj,t) - z x (xj,t). 

Note that in 0, B Xj ,Xj £ [0,1], represents the point-wise control located on the 
so boundary or in the domain. 

The space of weak solutions to System 0 is chosen to be fT 1 (0,1). Note that 
System 0 is exponentially stable in H [ (0.1) if the boundary controls Bq and B\ are 
chosen such that ko > 0 , k\ > 0, and ko + k\ >0 1231 . 

Denote a set of reference signals by {z[ > (xi,t)}^L 1 . Let 

e,(f) = z(xi,t) - zB{ Xil t) 

be the regulation errors. Let e(i) = {ei(t )}™ =1 and u(t) = {iij(t)}£L!■ 

55 Problem 1. The considered regulation problem is to find a dynamic control u(t) such 
that the regulation error satisfies e(t) —t 0 as t —> oo. 

The above in-domain control problem can also be formulated in another way by 
replacing the jump conditions in (|TT|) by point-wise controls as source terms. The 


resulting system will be of following the form 

m 

z t (x, t ) — z xx (x, t) = 5(x — 0 < x < 1, t > 0, (2a) 

j =i 

z(x, 0) = fi(x), (2b) 

Bqz = z x ( 0 , t) - k 0 z{ 0 , t) = 0, (2c) 

B\Z = z x (l, t) + k\z{l, t) = 0, (2d) 


where S(x — xfi) is the Dirac delta function supported at the point x 3 , denoting an 
actuation spot, and aij : t >->■ M, j = 1,..., m, are the in-domain control signals. 
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Lemma 1. Considering weak solutions in H 1 ( 0,1), System ([T]) and System 0 are 
equivalent if 


aj (t) = -Uj(t) = - [z x ] x . , j = 1, 


65 


PROOF. The proof follows the idea presented in (26). Indeed, it suffices to prove 
“System 0 =>■ System qr Let X = L 2 (0,1) be a Hilbert space equipped with the 
inner product ( u, v) = / u(x)v(x)dx. Let the operator A be defined by Au = u xx , 

Jo 

with domain T>(A) = {u £ H 2 ( 0,1); Bqu = B\u = 0}. It is easy to see that A*, the 
adjoint of A, is equal to A. Let A be an extension of A with domain 'D( A) = { u £ 


O '' b -L 

X ; u G U (xi,x i+ i)), B 0 u = Biu = 0, u(xj) = u(x~), j = 1, Let 

2—0 J J 

u G 'D(A), v G V(A*) = V{A). Using integration by parts we obtain that 


(Au,v) = ( u,Av) + ^2 ( u A x j ) - u *( x j~)) v ( x j)- (3) 

i=i 

Let X_i = ('D(A*)y, the dual space of T)(A). We need to define another extension 
for A. Let A : H 1 ( 0,1) —y X_i be defined by 


(Au, v) = (u, A*v) for all v G V(A*), 


(4) 


with V(A) = H 1 ( 0,1). Note that S(- — Xj) is not in X, but it is in the large space 
X_ -\. It follows from 0- 0. and A = A* that 

m 

Au = Au + ^2 i u x( x j) ~ u x (x+)) S(x - xj), ( 5 ) 

j=i 

in X _]. If u satisfies System 0^ then u(f) = Au(t), which yields, considering 0, 
ft(f) = Au + jyjL 1 (u x (x~) — u x (x+))5(x — Xj). Finally, we can see that System 0 
70 becomes System 0 with ay (t) = —Uj(t) = — [z x ] x . , j = 1, where we look 

for generalized solutions u(x , t) G T>(A) = H 1 ( 0,1) such that 0 is true in X_i. □ 


To establish in-domain control at every actuation point, we will proceed in the way 
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of parallel connection , i.e., for every Xj £ (0,1), consider the following two systems 


Zt(x,t) - z xx {x,t) = 0, x £ (0 ,xj) U (xj, 1), t > 0, 
z(x, 0) = <j>j(x), 

B 0 z = z x (0,t) - k o z(0,t) = 0, 

Biz = z x (l,t) + kiz(l,t) = 0, 
z(x+) = z(xj), 

B Xj z = [z x ] x . = Vj(t). 


(6d) 


(6b) 


(6a) 


(6c) 


(6e) 


(6f) 


and 


Z t (x,t) - z xx {x,t) = 5{x - Xj)fa(t), 0 < X < 1, t > 0, 
z(x, 0) = </>j(x), 

B 0 z = z x ( 0, t) - k 0 z{ 0, t) = 0, 

B 1 z = z x { 1, t) + k\z[l, t) = 0, 


(7d) 


(7b) 


(7a) 


(7c) 


with 4>j( x ) = <(>(x). Similarly, System (|6| and (|7]> are equivalent provided z £ 

H \ 0,1) and fa = -Vj = - [z x ] x .. Let aj = -Uj = fa = -Vj = ~[zi\ x . for any 
j = 1,2,to, where z 3 denotes the solution to System 0 One may directly check 
that z(x, t ) = i ^ ( x > t) ' s a solution to System Q. Moreover, 


m 



for any i = 1,2, Hence z(x,t) = YYjLi z-*(x,t) is a solution to System Q. 

Therefore, throughout this paper, we assume ctj = —Uj = fa = —Vj = — [z£\ x . for 
an Y J = 1, 2 ,to. Due to the equivalences of System 0 and Q, and System ([6]! and 
75 0 . we may consider Q and System ([6]) in the following parts. 

3. Control Design Based on Zero-Dynamics Inverse and Differential Flatness 

In the framework of zero-dynamics inverse, the in-domain control is derived from 
the so-called forced zero-dynamics that are constructed from the original system dy¬ 
namics by replacing the control constraints by the regulation constraints. To work with 
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the parallel connected system 0- we first split the reference signal as: 


z D {x,t) = ^2'y j (x,x j )zf(x j ,t), ( 8 ) 

i=i 

in which the function 7 j(x,Xj) will be determined in Proposition [4] (see Section |4j). 
Denoting by e J '(f) = (xj, t) — Zj(xj,t) the regulation error corresponding to Sys¬ 
tem 0 and replacing the control constraint by £ :l (t) = 0, we obtain the corresponding 


zero-dynamics for a fixed j: 

= £xx(x,t), x e (0 ,Xj) U (xj, 1), t > 0, (9a) 

£(*,0)=0, (9b) 

&(0,f)-fco£(0,i) = 0, (9c) 

^x(l) t) + &h£(l, t) = 0, (9d) 

€(xj,t) = Zj(xj,t). (9e) 


For simplicity, we denote by z ] and ^ the solutions to the j th systems 0 and 0, 
respectively. Also, we write henceforth Zj(t) = Zj(xj,t) as the reference trajectory in 
the j th system 0 if there is no ambiguity. The in-domain control signal of System 0 
can then be computed by 

v j = [ z x]xj = (10) 


Remark 1. Note that for ko > 0, k\ > 0 with ko + fci / 0, arguing as 4221/ . we have 
e J (f) = z^Xjjt) — Zj(xj,t) — »• 0 as t —» 00. 

Obviously, to find the control signals, we need to solve the corresponding zero- 
dynamics 0. For this purpose, we leverage the technique of flat systems EHED 
im In particular, we apply a standard procedure of Laplace transform-based method 
to find the solution to 0. Henceforth, we denote by f(x,s) the Laplace transform 
of a function f(x,t) with respect to time variable. Then, for fixed Xj £ (0,1), the 
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transformed equations of 0 in the Laplace domain read as 

s£(x, s) = lxx(x,s), x £ (0,Xj) U (Xj, 1), s £ C, (11a) 

OM) = o, (lib) 

f*(0,a)-fcof(0,a)=0, (11c) 

&(l,a) + fcif(l,8) = 0 ) (lid) 

£(xj,s) = Zj(xj,s). (lie) 

We divide o into two sub-systems, i.e., for fixed Xj £ (0,1), considering 

s£(:r, s) = £ xx (x, s), 0 < a; < Xj, s £ C, (12a) 

£(x,0)=0, (12b) 

fx(0,a)-fc 0 f(0,a)=0, (12c) 

€(xj,s) =Zj(xj,s), (12d) 

and 

s£(x,s) = £xx(x,s), Xj < x < 1, s £ C, (13a) 

£(x,0)=0, (13b) 

6e(1,s) + fci£(x,s) = 0, (13c) 

€(xj,s) = Zj(xj,s), (13d) 


Let £ J _(x, s) and £+(x, s) be the general solutions to ( | 1 2\ and ( p~3j ), respectively, 
and denote their inverse Laplace transforms by £?_{x,t) and £ J + (x,t). The solution to 
0 can be written as 

£ J (M) = €-{x,t)X{(p, Xj )} +C+(x,t)X{[x j ,l)}, 

where 

/ 1, x £ Qj C (0,1); 

I 0, otherwise. 

Then at each point Xi G (0,1), by © and the argument of “parallel connection ” (see 
Section^, we have [z x ] Xi = = [4]x ; = [4]x i5 i = Hence the 
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in-domain control signals of System o can be computed by 

^ 1, . . . ,777.. (14) 

In the following steps, we present the computation of the solution to System 
. Issues related to the reference trajectory z D (x, t ) for System G will be addressed 
so in Section [4] 

Note that £L{x, s) and £° + (x, s), the general solutions to ( fl2| and ( fT3| , are given by 

£t(x,s) = s) + C 2 <j> 2 {x,s), 

€+(x,s) = C 3 (/>i (x, s) + C4<j)2(x,s), 


with 


, , , sinh(v / sx) 

(pi{x,s) = - -= - 

Vs 

We obtain by applying ( |12c[ > and ( |12d[ ) 

Cl(j>l{Xj,s) + C 2 4>2(Xj,s) 

which can be written as 

IMxj,s) (j) 2 {xj,sy 
\ 1 -ko ) 


4 > 2 { x , s ) = cosh(i/sat). 


Zj(xj,s), Ci - k 0 C 2 = 0, 


^ = ^zf(x r: s) 


Let 


and 


R J _ 


( 4*1 ( x j > s) 

V 1 


02 (Xj,s) 
-ko 


Zj(xj,s) = — det(iV_)yl(a;.7, s). 

We obtain 

/C|\ _ adj(iV ) / zf(x. n s)\ _ I 

\Ca) det(f?i) \ 0 J y yL{xj,s) J 

Therefore, the solution to ( fl2| ) can be expressed as 

€-(x,s) = (k 0 (j>i(x) + ^ 2 (x))y°_(xj,s). 


(15) 


(16) 
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We may proceed in the same way to deal with (fi~ 3 j). Indeed, letting 


K = 


4 >i(Xj,s) (j) 2 (Xj,s) 

> 2(1 ,s) + fcl^l(l,s) s<j>i{l,s) + fci</) 2 (l, s) 


and 


z^(xj,s) = det (R J + )y J + (x J: s), 


( 17 ) 


85 we 


get from ([ 13 } 



(s<MM) +kihi 1 >s))y+(xj,s) 
-(</> 2 (l, s) + fci 0 i(l, s))yUxj, s) 


and 


\+{x, s) = ((s 0 i(l, s) + ki(/) 2 (l, s))</>!(x) 


+ (</> 2 (l,s) + k^^l, s))</> 2 (x)) y-[(xj, s). ( 18 ) 

Applying modulus theory |j 28 ll 29 l to ( [ 15 } and ( fT 7 } , we may choose yV x j , s) as the 
basic output such that 


y\{xj,s) = —det(R J _)yj(xj, s ), 

l+)yj( x j, 


V^-{xj,s) = det(Ri_)yj(xj, s). 


(19) 

( 20 ) 


Then, using the property of hyperbolic functions, we obtain from ( [ 16 } and ( [ 18 } that 

sinh(y / sa;j — y/s) 


¥-(x,s) = fci 


>/* 


— COSh(v / StCj — yfs) 


. , sinhf-v/stc) , , r- \ \ ~ 

x [ k 0 - -J= -b cosh(v sx) ) yj(xj,s), ( 21 ) 

co , y 1 , sinh (y/sx - y/s) 

£+0,s) = ( fci-- 

x ( k 0 ) +cosh(v / sa;j) ) yj(xj,s). ( 22 ) 


— cosh(v / sa: — V~s) 


Note that 


€ J (x,s) =V-(x,s)x{( o,x 3 )} + £+ 0 E > s )X{[x. h i)} 


( 23 ) 
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is a 


solution to (fIT|). Using the fact 


00 x 2n +i 


sinh.T = 7 -—, cosh a: = 

(2n + 1 )! Z^ 


n —0 


n —0 


( 2 n)! 


we obtain 


? (a:, s) 


' 00 n 2k+It 1 \2(n-fc) + l 00 « „2fc+l/'„ 1 \2(n-fc) 

koki EE - E... 


^ (2*+ 1)![2(»-*) + !]! 


n=ofc^o ( 2 & + l)![ 2 (n — &)]! 


*( x i -l) 2(n ~ k)+ \ n ST'ST' x 2k (xj - l) 2( "- fc) _ f 
+ fcl 2 ^ i9Mb9i'?7 _ m 4- ill s _ Z^ Z^ Z^\iTo 7 ZT 7 Tvrr s |X{(o,a^)} 


00 n 2ft c 


(2*)![2(n - k) + 1]! (2fc)![2(n - At)]! 


n—0 k—0 

oo n 2fc+l 
X 3 


fc o fc i E E 


atj" ‘ ~(x - i) 2 («-fe )+ 1 


00 n 2fe+l 
X 3 


(2k + l)![ 2 (n — k) + 1 ]! 


^»EE 


(a: - l) 2 ("- fc ) 


n—0 /c=0 

oo n „2k( x _ -Q2(n-fc)+l 


n—0k =0 


(2k + l)![ 2 (n — fc)]! 


+ ' tl SS(i)W-‘) + l]l S " SS ( 2 *) 1 [ 2 (»-k)]l 


oo n 2ft 
X 0 


(x- l)^ n ~ k ) 


s Xffo.i)} 


Vo¬ 


lt follows that 


f (M) 


, 0 1 SS ( 2 k + l)![ 2 (n - k) + 1 ]!^ 0 EE ( 2k + l)![ 2 (n - k)]\ V > 


^ UL 2,2ft 


+^EE 


(Xj - l) 2 ("- fc ) + ! ( . 


n) 


00 n ft 


n= 0 ft=o ( 2 A) ! [ 2 (ft — A) + !]• ^ ( 2 *)![ 2 (n-fe)]! 


-EE 


( x . _ l)2(n-ft) 


n) \ 

Vj X{(o,* 3 )} 


oo n 2/c+l 
X 3 


fc ofci 

V n—0 /c—0 


(x _l)2(n-ft) + l (n) ; ^^xf+^at-l^-ft) _ (n 

Vj '°EE ( 2 k + 1)![2 (n - k)]\ Vj 


(2k + l)![ 2 (n — fc) + 1 ]!' 


n—0 /c—0 


oo n 2fcf 
X 3 


1 SS (2fc) ! [2(«- As) + i]!^ EE ( 2 A)![ 2 (n-fc)]!^ J X{fe ’ 1)} _ 

(24) 


By a direct computation we get 

k 0 k ± 


P 

^>X 


+ y/s ) sinh(v / s) + (fc 0 + fci) cosh(y / i) 


s r ‘ 


x Vj(xj, s) 


°° s ra+l \ 

= (E pSTl)! + + *'>§« + S(STTj!) 


n=0 


( 2 n)! ' ( 2 n + 1 )! 

n—0 v J n—0 v ' ■ 


(25) 


11 




































It follows from (jT4j) that 


Uj(t) = [££], 

=k 0 ki ^2 


(r 

y) 


\xj,t) 


4 (2n + l)! 

n—0 7 


(/c 0 +fci) ^2 


n—0 


yj n) (xj,t ) 
(2n)! 


00 


E 

n—0 


(2n+ 1)! 


(26) 


Remark 2. 77ze reference signal Zj(xj,t) can be derived in the same way from the 
flat output from CD and O- However, as the flatness-based control is driven by flat 
output, there is no need to explicitly compute Zj(Xj, t). 


so 4. Motion Planning 

For control purpose, we have to choose appreciate reference trajectories, or equiva¬ 
lently the basic outputs. In the present work, we consider the set-point control problem, 
i.e. to steer the heat distribution to a desired steady-state profile, denoted by z D {x). 
Without loss of generality, we consider a set of basic outputs of the form: 

Vj{t) = y(xj)(pj(t), j = 1,. ..,m, (27) 

where <pj(t) is a smooth function evolving from 0 to 1. Motion planning amounts 
then to deriving y{xf) from z D (x) and to determining appropriate functions Pj{t), for 
j = 1,... ,m. 

To this aim and due to the equivalence of the systems ([TJ) and ([2]), we consider the 
steady-state heat equation corresponding to System 

m 

z xx (x) = ^ 2 S( x ~ x :i) (x :ii 0 < x < 1 , t > 0, (28a) 

i=i 

^(0)-M(0)=0, (28b) 

z x (l) + k 1 z{l) = 0. (28c) 

Based on the principle of superposition for linear systems, the solution to the 
steady-state heat equation ( |28| can be expressed as: 

m m 

z( x ) = ^2 g ( x ^O s (C - x j)oij d C = '52G(x,x j )a j . (29) 

J 0 j =1 j —1 
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where G(x. (j) is the Green’s function corresponding to ( |28j ), which is of the form 

(fctC - h - 1 )(k 0 x + 1) 

G(x, 0= ' 


kg + /Ci + 

(ftlX - fel - l)(fcoC + 1) 


0 < x < C; 


( 30 ) 


c < X < 1. 


kg + k\ + kgk\ 

95 Indeed, it is easy to check that G xx (x, () = S(x — () and G(x, () satisfies the bound¬ 
ary conditions, G^O, 0 — fcoG(0, £) = 0 and G x (l,£) + &iG(l,£) = 0, the joint 
condition, G(( + ,() = an d the jump condition, [G x (x, C)]f = 1. 

Taking m distinguished points along the solution to ( [28] ), z(x\,t ),..., z(x m , t), 
we get 

^ z(xi) \ f G(x i,xi) ••• G(xi,x m )\ f a 

; . od 

\z(Xrn) J \G{Xjn^X\) *** *£ra)/ \&m J 


Lemma 2. The matrix (G(xi, Xj)) mxrn chosen as in © is invertible. Thus 


(ar\ 


(G(x 1 ,x 1 ) ■■■ G(xi,x m ) \ (z(x i)\ 


(32) 


\pim J yG(^r m ,x i) ••* G(x m ,x m )y 'y*'{x rn ) J 

Proof. For m = 1, since kg > 0, k\ > 0, fco + k\ > 0 and x\ € (0,1), it follows 
that G(xi, x±) = fc 1 x ^ 1+1 ^ < 0- Hence it is invertible. We prove the claim 
for to > 1 by contradiction. Suppose that the matrix ( G(x-i, Xj)) mX m is not invertible, 
then it is of rank to — 1 or less. Without loss of generality, we may assume that for 
some x n > Xj with i = 1 ,..., n — 1, there exist n — 1 constants Z l5 Z 2 , l n -1 such 
that 


n—1 


G(xi,x n ) = hG(x 1 ,x i ), 

i= 1 


(33a) 


G{x n ^x n ') — ^ ^ (zG(x n , Xj), 


(33b) 
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where 1 < n < m and Y^i=i > 0- Let 


G(x) = G(x,x n ), F(x) = ^ hG(x,Xj). 


i =1 


ioo ([33} shows that F(x) = G(x) at every boundary point of [ 21 , 22 ], [ 22 , £ 3 ], ■ • ■, 

Note that F(x) is a linear function in [ 21 , 22], [22,23], , [ 2 „_i, x n ], and that G(x) = 
^ X ko+kr+lo l-’T 1 ~ [^ti^n]- i- e -’ G(a;) is a linear function in [ 21 , x n ]. Hence 
F( 2 ) = G{x) in [ 21 , 2 „]. 

By F( 2 i) = G( 2 i), we get 

n— 1 

k±x n - fci - 1 = li(kiXi - hi - 1). (34) 


By F{Xji) = G( 2 „), we get 


j=l 


n—1 


^ 0 "T 1 ^ ^ + l). 


Therefore 


n— 1 


X> = 1 - 

By F x ( 2 f) = G x ( 2 f) and F x ( 2 “) = G x ( 2 “), we get 

n— 1 

k 0 {kix n - kx - 1) =fc 0 /ci liXi - k 0 (ki + 1) li 4- Zifci, 


(35) 


(36) 


n—1 


2=1 
n—1 


n—1 


=k 0 ki yy Zj2j+fci y^ k. 


(37) 


It follows that h = 0’ which yields, considering ([36}, l\ = 1. By F x ( x 2 ) = 

Gx ( 22 ) = F x (x 2 ), we deduce 


liki(k 0 xi + 1) + hki(koX2 + 1) + k 0 y^ U(k\Xi — k\ — 1) 

i=3 

n— 1 

=hki(k 0 xi + 1) + fc 0 ^ k{k\Xi -k\- 1), 


i =2 


which gives l 2 = 0. Similarly, by F x (2+) = G x (xj) = F x (x j ) (j = 3,4,..., n - 2), 
we obtain (3 = (4 = ... = l n -2 = 0. Hence Z n _i = 0. Then we deduce from ( |37} that 

k 0 {kix n — fci — 1) = fci(fe 0 2 i + 1). (38) 
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It follows from ( |35[ > that 


fco^i + 1 = k 0 x n + 1. (39) 

We conclude by ( |38| ) and ([39]) that k 0 + k-\ + k 0 k\ = 0, which is a contradiction to 
105 ko >0, k\> 0, and fco + k\ > 0 . □ 

In steady-state, we can obtain from ( [26] ) that 

Uj = (koh + k 0 + ki)y(xj) = -otj. (40) 

Finally, y(xj) can be computed by ( [32] ) and ( [40] ) for a given z D (x). 

It is worth noting that ( [32] ) provides a simple and straightforward way to compute 
the static control from the prescribed steady-state profile. Indeed, a direct computation 
no can show that applying ( [32] ) will result in the same static control obtained in ll22ll where 
a serially connected model is used. 

To ensure the convergence of ( |24| and ( |26[ ). we choose the following smooth func¬ 
tion as (pj{t): 

if t < 0 

if t e (0, T) (41) 

if t > T 

which is known as Gevrey function of order rr - 1 + 1/c, c > 0 (see, e.g.. Ill 3fl ). 


0 , 


= < 


[ exp(-l/(r(l - r))) E dr 
Jo _ 

[ exp(—1/(r(l — r))) e dr 
Jo 


Lemma 3. If the basic outputs j = 1,..., m, are chosen as Gevrey functions of 

order 1 < o < 2, then the infinite series m and (]26[) are convergent. 


PROOF. We prove the convergence of the power series ( |24| ) and ( [26] ) using Cauchy- 
Hadamard Theorem. Indeed, it suffices to prove the convergence of 




(42) 


ns Then the convergence of the series in ( [24] ) and ( |26| ) follows easily using the same 
argument. 
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We recall that the bounds of Gevrey functions of order a are given by 


3K,M > 0,Vfc G Z> 0 ,Vi g [f 0 ,T], 




(fe+i) 


(t) 


< M 


(vr 

I\ k 


( 43 ) 


Denote in 


bn=J2 


1 


k =0 


(2fc)![2(n - k)\\ 3 


vT'm- 


Then, ( |42| converges if limsup^^^ < 1. Now b n can be estimated by ( |43f > 

M [n\y 


I K\ < E 


k=0 


< M 


(2k)\[2(n - k)\] K n 

2 n ( n\y 


(2n)! AT" ' 


Therefore 


__ 9 , f (n l ) 1 / rl ') Cr 

limsup l/EI < limsup — M 1 ^— K — - 

n —roo n—too (((2d)!) 

2 (n/e) CT 


< limsup 


K (2n/e) 


2A-““ S T"" = i 2A’ ” 2 ’ 

cr > 2, 


0, 

1 

2iT ’ 
oo, 


(T < 2, 


where in the second inequality we applied Stirling’s formaula yn! ~ {n/e). We can 
conclude by Cauchy-Hadamard Theorem that ( |42[ ) converges for a < 2, and for cr = 2 
if 2K > 1. The series ( |42j i diverges if cr > 2. □ 

Theorem 4. Assume fc 0 > 0, fci > 0, with k 0 + k\ > 0. Let the basic outputs 
j = 1,... ,m, be chosen as ED with an order 1 < cr < 2. Let the reference trajectory 
of System 0 be given by 0 with 


7 j(x,Xj) = - 


(fcofci + k 0 + ki)G{x,Xj) 


, j = 


(44) 


{k 0 Xj + 1 ){k 0 (xj - 1) - 1) : 
where G{x,(j) is the Green’s function defined by Then the regulation error of 

System 0 with the control given in ( |26| l tends to zero, i.e., 

efit) = z{xi, t) — zf (xi, t) —>■ 0 as t —)• oo, 


120 for i = 1, 2, m. 
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PROOF. By a direct computation we have 


\ei(t)\ =\z(xi,t) - z^(xi,t)\ 


(k 0 k! + k 0 + k 1 )G(x i ,Xj)z^{xj,t) 
+ Z^ (k oXj + l)(k 0 (xj - l) - 1) 

(k 0 k i +k 0 + ki)G(xi, Xj)£ J {xj,t) 


.7 = 1 


z(xi,t) + ^2 


(k 0 Xj + l)(k 0 {xj - 1) - 1) 


< | z(xi,t) - z(xi )| + 


+ (fcofct + k o + fcr) ^ G(xj, Xj)y{xj) 


7=1 


/ 7 7 ,7 , 7 X V' G^X^^X^t) 

0 ,1 ° 1 ^ (fco^-+ l)(fco(a;j - 1) - 1) 


— (AsoAsi + fc 0 + k 1 )y^G(x i ,Xj)y{xj) 

7 = 1 


By ( [40] ) and ([29]), it follows 


z(xi) = ~(k 0 ki + fc 0 + k\)' s ^^G{x i ,Xj)y{xj). 

7=1 

Based on {24], ( |27] >. and the property of ipj (t) we have 

G(xj,t) _ 

7 -- w , . -—-r —>■ via;,) as f —>• 00 . 

(k 0 Xj + 1 )(k 0 (xj - 1 ) - 1 ) 

Note that z(xi, t) —> z{xf) as t — > 00 . Therefore \ei{t)\ — > 0 as t — t 00 . 


□ 


Remark 3. For any x £ (0,1), replace x, by x in the proof of Proposition^ we can 
get \ z(x, t) —z D (x, t)\ —> 0 as t —> 00 , which shows that the solution z(x, t ) of System 
Q converges to the reference trajectory z D {x, t) at every point x £ (0,1). 


125 5. Simulation Study 

In the simulation, we implement System ([2]) with 12 actuators evenly distributed in 
the domain at the spot points {1/13, 2/13,..., 12/13}. The heat distribution and the 
time are all represented in normalized coordinates. The numerical implementation is 
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based on a PDE solver, pdepe, in Matlab PDE Toolbox. 200 points in space and 50 
130 points in time are used for the region [0, 1 ] x [0, 2] in numerical simulation. The basic 
outputs ifj (t ) used in the simulation are Gevrey functions of the same order. In order to 
meet the convergence condition given in Proposition [3j the order of Gevrey functions 
is set to a = 1.1. The feedback boundary control gains are chosen as = k\ = 10. 
The initial condition in simulation is set to z(x, 0) = cos(7rai). 

135 The desired heat distribution and the evolution of heat distribution of the system 
controlled by the developed algorithm are depicted in Fig. [T] Snapshots of regulation 
errors are presented in Fig. [2] which show that the regulation error tends to 0 along the 
space. The control signals that steer the heat distribute from the initial profile to track 
the de sired one are illustrated in Fig. [3] The simulation results show that the system 
uo performs very well with affordable control efforts. 



(a) 


(b) 


Figure 1: Evolution of heat distribution: (a) desired profile; (b) solution surface. 




(a) (b) 


Figure 2: Snapshot of regulation errors: (a) errors for t £ [0,1]; (b) errors for t £ [1,2]. 
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(a) 


(b) 


Figure 3: Control signals: (a) u\ to Mg; (b) M 7 to M 12 . 


6. Conclusion 

This paper presented a solution to the problem of set-point control of heat distribu¬ 
tion with in-domain actuation described by an inhomogeneous parabolic PDE. To apply 
the paradigm of zero dynamic inverse, the system is presented in an equivalent parallel 
us connection form. The technique of flat systems is employed in the design of dynamic 
control and motion planning. As the control with multiple in-domain actuators results 
in a MIMO problem, a Green’s function-based reference trajectory decomposition is 
introduced, which considerably simplifies the design and the implementation of the 
developed control scheme. The convergence and solvability analysis confirms the va- 
150 lidity of the control algorithm and the simulation results demonstrate the viability of 
the proposed approach. Finally, as both ZDI design and flatness-based control can be 
carried out in a systematic manner, we can expect that the approach developed in this 
work may be applicable to a broad class of distributed parameter systems. 
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